Dynamical Franz-Keldysh effect 
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We introduce and analyze the properties of dynamical Franz-Keldysh effect, i.e. the change of 
density-of-states, or absorption spectra, of semiconductors under the influence of time-dependent 
electric fields. In the case of a harmonic time-dependence, we predict the occurence of significant 
fine structure, both below and above the zero-field band-gap, which should be experimentally ob- 
servable. 
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Almost forty years ago Franz and Keldysh |] 
pointed out that static electric fields modify the linear 
optical properties of bulk semiconductors near the opti- 
cal absorption edge: the absorption coefficient becomes 
finite (though small) for photon energies below the band 
gap, and that the absorption coefficient shows oscilla- 
tions for energies above the band gap. Studies of this 
kind regained importance with the development of quan- 
tum confined semiconductor structures, and a large body 
of both experimental and theoretical knowledge has been 
developed since mid 80 's ||. The purpose of this Letter 
is to point out that the absorption curve develops addi- 
tional structure when the sample is placed in a strong 
ac-field, such as generated by a free-electron laser. As we 
shall demonstrate below, the zero-field density of states 
splits into replicas, with field-dependent shifts and time- 
dependent amplitudes. In contrary to the static case, 
where it is sufficient to solve a simple Schrodinger equa- 
tion, the time-dependent case necessitates the use of a 
more advanced formalism; in the present case we apply 
nonequilibrium Green function techniques. In order to 
present our predictions in a clear-cut manner, we have 
stripped the model calculation from complications that 
may be present in real samples, such as scattering due 
to sample nonidealities and phonons, or excitonic effects. 
Thus, the main motivation of the present analytic study 
is to point out the existence of an effect, and not address 
the quantitative details, which can be treated properly 
only with a full-scale numerical calculation based on, for 
example, solution of the semiconductor Bloch equations 

§• 

The calculation of the electroabsorption coefficient 
a(u>) for free carriers is discussed in text-books ||, and 
the result can be expressed as 



(1) 



where vb n and E n are the eigenfunctions and eigenvalues 
of the electron-hole pair Schrodinger equation, respec- 
tively. The important feature of Eq. ([j]) for our purposes 
is that it states that the absorption coefficient is essen- 
tially given by the joint density-of-states, a{w) ~ p(ui) = 



^ n 5{u> — E n ). In equilibrium many-body theory one 
can express the density-of-states in terms of the imagi- 
nary part of the retarded Green function, or, equivalently, 
by the spectral function. The important property of the 
nonequilibrium theory of Kadanoff-Baym || and Keldysh 
is that in nonequilibrium the diagrammatic structure 
of the perturbation theory is the same as in equilibrium, 
with the proviso that real-time integrals must be gen- 
eralized to a complex contour j^]. The upshot is that 
the nonequilibrium retarded Green function can be cal- 
culated from a Dyson equation which is formally identical 
to the equilibrium case. Thus, we can analyze the elec- 
troabsorption properties of a semiconductor in an applied 
time-dependent electric field by evaluating 



p(u,T) = — V ImG r (k,w,T) 

k 

= ^-£I(k,a,,T). 



(2) 



Eq.(|2|) introduces a generalized density-of-states function, 
which is time-dependent due to the external electric field. 
Whether a measurement probes p(uj, T) or perhaps its 
suitably weighted time-average depends on the various 
time-scales of the particular experimental set-up; here 
we focus on the full time-dependence, because it con- 
tains most information. In (||) we also introduce the 
gauge-invariant form G r (k, u>,T) of the two-time Green 
function G r (x, x', t, t'), which we calculate via its Dyson 
equation. Effects due to interactions can be taken into 
account by introducing a suitable self-energy. Under 
certain conditions the Dyson equation for the retarded 
Green function may involve particle densities; in this 
case one must also solve the appropriate quantum ki- 
netic equation for the distribution function. However, 
this problem is not addressed in the present paper. 

In an actual calculation the longitudinal electric field 
must be introduced by using some explicit gauge; how- 
ever predictions for observables must clearly be gauge- 
independent, and hence the gauge-invariant formulation 
is essential. In general, the transformation to gauge in- 
variant variables occurs via 191 
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FIG. 1. The function Ai as a function of the variable 
uj* — (lu — £fc)/f2. The field strength is given by luf/SI = 1. 
The 5-function singularities of A\ have been broadened for 
illustrative purposes 



G r ( P ,X)= / AG r (g,X) 



exp(ig M [p p 



1/2 



1/2 



dXA^X + Xq)}) , (3) 



where we have used the Einstein summation conven- 
tion and a "covariant" notation, where four-momenta 
and coordinates are given by p = (w,k),<7 = (r, r) = 
(t-t',x-x'), and X = (T,R) = ([t + £'])/2, [x + x']/2). 
Correspondingly, the potentials defining the electric field 
are contained in A = (</>, A). 

In our calculation of the density-of-states for a uni- 
form but time-dependent field it is convenient to de- 
scribe the electric field with a vector potential, A(t) = 
— E s'm(nt)/n (here we restrict ourselves to harmonic 
time-dependence). The retarded Green function satisfies 



- e[p - A(t)] }G r (p, i, O = *(t - f) 



(4) 



Using (||) Eq. (|]) immediately leads to the gauge invariant 
spectral function A = i(G r — G a ), 

A(k,u,T) = J drdre tw J dpe lr p 

(•T+t/2 
IT-t/2 

where 



r r T+T/2 ^ 

expj-z/ dhe[p- A(h)] \, (5) 

1 JT-t/2 ' > 



to = tw — r ■ k — / — r-A(fi) 

JT-t/2 t 



(6) 



It is illustrative to see how the static field results |Tl]] 
can be recovered from Eqs.(|^-||). In the limit £1 — > and 
for parabolic bands Eq.rta) reduces to 



FIG. 2. The function A2. The parameters are defined in 
the caption of Fig. 1. 



A dc (k, t) = ex P H(e fe T + E 2 t 3 /2A)} 
A dc (k,u) = (2n/0)Ai[{e k -u)/0\ , 



(7) 



with = (h 2 e 2 E 2 /8m) 1/3 1 and Ai(x) is the Airy func- 
tion. Evaluation of the momentum sum leads then in 
three dimensions to 

P %(u>) cx [Ai' 2 (-(^/e)) + (c/e)Ai 2 (-(c/e))] , (8) 

with 9 = 4 1 / 3 /3/?i. Eq.(g) is the standard form for the 
Franz-Keldysh line-shape ||. 

In the time-dependent case, one finds 

A(k,cj,T) = J dTCOs[(uj-e k -LU F )T + X(T) + 2Y(T)] , 

(9) 



where we have introduced the notation 
X(t) = -jj- sin Qt cos 2QT 



Y(r) 



ujf 4sin 2 (fir/2)sin 2 fiT 



fir 



(10) 



(11) 



and defined the field parameter up = e 2 F 2 /(4mM7 2 ) 
|]l2f . An explicit analytic evaluation of the integral in (|]) 
is not possible, nor is a straightforward numerical cal- 
culation feasible due to its singular character. The sin- 
gularities, however, can be isolated with analytic means, 
and we proceed as follows. The trigonometric identity 
cos(a; + y) = cosx — 2sin(ir. + y/2) sin(y/2) allows us to 
write (||) as a sum of two terms, A = A\ + A2, where 



A 1 (k,w,T) = J dTCOs[(Lo-e k -u F )T + X{T)] , (12) 

^(k, lu, T) = — 2 J dr sin [(a; — — lof)t 

+X(t) +Y(t)] sinF(r) . (13) 




FIG. 3. Density of states, three dimensions. 

The integral defining A 2 is convergent, and its numer- 
ical evaluation is straightforward. On the other hand, 
analytic progress can be made with A±, which contains 
the aforementioned singularities. Recalling the identity 
f^° oo dxcos(ax - bsxnx) = 2wY,™ = - 00 J n (b)S(n - a) , 
where J n is the n-th order Bessel function, we finally 
obtain 

oo 

Ai(k,w,T) = 2»r Jn[-^cos(2fiT)] 

n— — oo 

x 5(u> - e k —nVt — uj f ) . (14) 

These two components A\ and A2 of the spec- 
tral function have quite distinct properties. The fre- 
quency sum rule is exhausted by A\ alone, jl3| i.e. 
J dw/2wAi(uj,T) = 1, while / cIuj/2ttA 2 (lu, T) = 0. Fur- 
ther, the delta-function singularities of A\ are reminis- 
cent of 'photonic side-bands JL4| , however in the present 
case these features are shifted by the field- and frequency- 
dependent parameter u)p. In the zero-field limit A\ re- 
duces to the field-free result ao(k, uj) = 2nS(uj — ek), while 
A 2 vanishes. Note also that the spectral function de- 
pends quadratically on the electric field. This is consis- 
tent with physical intuition: since there is no preferred 
direction, it should not matter if the field direction is re- 
versed, E — > — E, and thus one expects to find a spectral 
function that is even in the applied field. The functions 
A\ and A 2 are displayed in Figs. 1 and 2, respectively. 

We can use the above results to calculate a time- 
dependent density-of-states function p(uj,T). The con- 
tribution arising from A\ is readily evaluated by using 
the delta-function, while A 2 requires some more work. 
In three dimensions we find |l5|] 

2ir 2 h „ = _ DO si 
x (u-uj f -nsl) 1/2 (15) 
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FIG. 4. Density of states, two dimensions. 

and 

x sin[(w - uj f )t + X(t) + Y(t) + vr/4] . (16) 

Here the summation over n should be cut-off at the 
largest integer not exceeding [(uj — ujf)/Q]. A numeri- 
cal example is shown in Fig. 3. The overall qualitative 
picture is dominated by the free-electron like square-root 
components generated by the photonic side-bands in A\ . 
At different points of time, however, the n:th side-band 
amplitude, which is determined by the Bessel function 
J n , may change its sign, which can be seen in Fig. 3. 
The more detailed fine structure in p 3D (uj,T) is caused 
by the contribution originating from A 2 . 

Similar, but even more dramatic effects can be found 
in two dimensions where the result is 

p^( w> D = -^j ± Jn [-^ cos (nT)] 

Znfi si 

n— — 00 

x 9(lu - u> F - nQ) (17) 

#>,T)= ™ r^ S m[Y(r)] 
2n l h J-oo t 

x cos[(uj-uj f )t + X(t) + Y(t)} . (18) 

The zero-field step-function reflecting the two dimen- 
sional density-of-states is split into several substeps ex- 
tending to energies well below the gap, and the above 
threshold values are modified by additional structures, 
whose amplitude tends to zero for increasing energies. 
Just as in the three-dimensional case, the relative ampli- 
tude of the substeps changes as a function of time, and 
the fine structure in the energy dependence is due to A 2 . 

It is also of interest to examine how the time-dependent 
density of states depends on the field strength at a 
fixed frequency. The effects are most easily analyzed 
by examining the time- average or the density of states, 
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FIG. 5. Time- averaged 2-dimensional density of states (in 
units of zero-field density of states) for ujf/Q.= 0.5 (dashed 
line), 1 (solid line), and 2 (dash-dots). 

p avo (w) = (0/2tt) J^ /Q dTp{u,T), and Fig. 5 shows the 
results of such a calculation. The following features of 
p ave (uj) are of importance: (i) The main absorption edge 
moves to higher energies with increasing field strength; 
(ii) There is a significant enhancement at the absorption 
edge; and (iii) the range of energies where the edge is 
modified increases with increasing field strengths. 

Finally let us comment on the observability of the pre- 
dicted effects. There are several limiting factors: the 
spectral resolution in an absorption measurement, the 
available frequencies, and the maximal ac-field strengths 
that can be applied to the sample. From the results pre- 
sented above we conclude that most favorable conditions 
are obtained if cjp/f2 > 1. It would appear that this con- 
dition is most easily achieved at low frequencies. How- 
ever, a small Q implies that all the fine-structure would 
be confined to the immediate neighborhood of the ab- 
sorption edge (in the example of Fig. 5 the fine structure 
extends only a few O's above the edge), and not resolv- 
able spectroscopically. Increasing fi leads to an increased 
E in order to maintain sufficiently large luf/Q. By us- 
ing parameters attainable with free-electron lasers |fl6| , 
Q ~ 1 THz, E ~ 1 MV/m, and the GaAs effective mass, 
we estimate that ujf/^I ~ 1. This implies that the fine- 
structure extends a few meV around the absorption edge, 
and should be observable. 

In summary, we have calculated the density of states of 
free carriers in time-dependent fields, and find that a har- 
monically varying external field leads to a field and fre- 
quency dependent shift of the main absorption edge, and 
gives rise to fine structure ("dynamical Franz-Keldysh 
effect") of observable magnitude. 

We acknowledge useful comments from Prof. S. Koch 
and Prof. Jai Singh. 
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